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Abstract 

We study the neutralino relic density in the minimal supersym- 
metric standard model with natural light Higgs sector in which all 
r*j . Higgs masses, the supersymmetry (SUSY) breaking parameters, and 

Q-f the higgsino mass parameter \i are of order the weak scale. To realize 

O ,■ this situation we adopt nonuniversal Higgs masses at the grand uni- 

fied scale. We show that in some parameter space in which the SUSY 
breaking parameters are comparatively small, not only the constraint 
from the observed relic density of dark matter but also the LEP Higgs 
►> . bound and the constraint from the b — > 57 process are satisfied. In 

most of the parameter space, the neutralino relic density becomes 
smaller than the observed relic density in contrast with the results in 
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(f-) ■ the constrained minimal SUSY standard model (CMSSM). The reason 

_il , is that the neutralino coannihilation processes to Higgs bosons open 

^^ ' even if the gaugino mass is small and the cross sections become large 
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due to the small dimensionful parameters. Especially small /i param- 
eter and the light CP-odd Higgs, which are difficult to be realized in 
the CMSSM, are essential for the result. 
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1 Introduction 

The minimal supersymmetric (SUSY) standard model (MSSM) is one of the 
hopeful extensions of the standard model (SM). It is attractive not only in 
the point of the weak scale stability, but also in the fact that SUSY models 
with R-parity have the lightest supersymmetric particle (LSP) as a good 
candidate for dark matter. Since the thermal relic density of the LSP can 
be calculated once we fix the parameters in the MSSM, it is interesting to 
examine parameter space, which is consistent with the observed value 

n DM h 2 = n m h. 2 -n b h. 2 (i) 

= (0.1277+°;°^°) - (0.02229 ± 0.00073) , (2) 

where $lDM,^m an d &b are the energy densities of dark matter, matter, and 
baryon of the universe [Tj. Here, h is the normalized Hubble parameter 
such as the present Hubble constant is given by Hq = 100hkms~ 1 Mpc _1 . 
Many studies have been done about the neutralino relic density in the con- 
strained MSSM (CMSSM) [2], in which all dimensionful parameters can be 
presented by only five parameters, the unified gaugino mass mi/2, the uni- 
fied scalar mass mo, the universal couplings for three scalar vertex A, the 
parameter for Higgs mixing B, and the higgsino mass /i. Unfortunately the 
allowed region for the parameters in the CMSSM are quite limited because 
in most of the parameter region consistent with experiments, the calculated 
thermal relic density of the neutralino become too large to satisfy the ob- 
served value. Moreover, in the CMSSM, the LEP constraint to the standard 
model Higgs mass, rrih > 114.4 GeV (95% C.L.) [3j, requires comparatively 
large SUSY breaking parameters in order to make the lightest MSSM Higgs 
heavier via loop corrections than the upper bound for the SM Higgs mass 
nih > 114.4 GeV. Such large SUSY breaking parameters destabilize the weak 
scale. This problem is called as the little hierarchy problem. 

Recently, it has been pointed out that in the nonuniversal Higgs mass 
model, the LEP constraints can be avoided due to the smaller ZZh coupling 
than in the SM [HElElEj- Here, h is the lightest CP even Higgs. Therefore, 
the large SUSY breaking parameters are not needed. This avoids the little 
hierarchy problem. To obtain the small ZZh coupling, generically, the light 
Higgs sector is required, in which not only the usual CP even Higgs but 
also the other Higgs bosons have the weak scale masses. Contributions of 
the light charged Higgs to the b — > 57 process would be too large to be 
consistent with the experimental value, if the chargino contribution has not 
compensated the charged Higgs contribution. This cancellation due to the 
supersymmetry works well [El El QUI E] because all the mass scales in the 



Feynman diagrams contributing the b — ► S7 process are of order the weak 
scale in the models with the natural SUSY breaking parameters and the light 
Higgs sector. In such models all the dimensionful parameters are of order 
the weak scale. The charged Higgs mass m#± and the Higgsino mass fi are 
fixed at the weak scale. The sfermion mass mo, the gaugino mass mi/2 and 
the scalar three point coupling A are fixed at the GUT scale. We call such 
scenario the natural light Higgs scenario. 

In this paper, we calculate the thermal relic density of the lightest neu- 
tralino in the natural light Higgs scenario. The result is totally different 
from the result in the CMSSM. The thermal relic density in this scenario 
tends to be smaller than the observed dark matter energy density. This is 
mainly because the neutralino coannihilation modes to Higgs bosons such as 
X1X1 — *■ hA and x?X? — > HA open due to the light Higgs sector and because 
the cross sections become large due to the small dimensionful parameters, 
especially small \x. Here Xi: A an d H are the lightest neutralino, the CP 
odd Higgs, and the heaviest CP even Higgs, respectively. The larger total 
annihilation cross section of the neutralino leads to the smaller thermal relic 
density. If // is large, the cross sections of neutralino coannihilation processes 
to Higgs bosons decrease because of the small Higgsino components of the 
lightest neutralino. Thus the energy density of the neutralino becomes larger 
than that with small fi. 

There are two essential points. One of them is that the light Higgs 
bosons make it possible to open the neutralino coannihilation processes to 
Higgs bosons. And the other is that the small \x parameter makes the 
cross sections large. In the CMSSM, it is difficult to satisfy both of them. 
Roughly speaking, to obtain the light Higgs sector, the both mass parameters 
mi 
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must be around the weak scale, which is difficult to be satisfied in the CMSSM 
because m 2 — m 2 , becomes much larger than the weak scale. Here, H u and 
H d are up-type Higgs and down-type Higgs, respectively. Actually, in the 
CMSSM, the difference m\ — m 2 at the weak scale becomes roughly Sm 2 , 2 , 
where rrii/2 is taken roughly larger than 300 GeV to satisfy the LEP constraint 
to the SM Higgs mass bound. 

For the reasons stated above, the neutralino relic density in the natural 
light Higgs scenario tends to be smaller than the observed energy density. 
And, in this scenario, there are parameter regions where the neutralino relic 
density agrees with the observation. 



In section [21 we show the numerical calculation of the neutralino thermal 
relic density in the natural light Higgs scenario. After a discussion about the 
allowed region, we conclude in section El 

2 Neutralino relic density in natural light Higgs 
scenario 

In this section, we calculate the neutralino thermal relic density numerically 
in the MSSM with the light Higgs sector and natural SUSY breaking param- 
eters. There are two additional dimensionful parameters, m 2 H and m 2 H in 
the nonuniversal Higgs mass model. Then we have seven parameters. One 
of seven parameters is fixed by the Z boson mass, and thus we have six pa- 
rameters. In this paper, three of them, the universal sfermion mass mo, the 
gaugino mass m^, and the universal coupling for the three scalar interaction 
Aq are fixed at the GUT scale, and the other parameters, /i parameter, the 
ratio of two Higgs vacuum expectation values tan /3, and the CP odd Higgs 
mass rriA are fixed at the weak scale. 

In this paper, we adopt the small ZZh coupling, scenario [U El El E] to 
satisfy the LEP constraints to the SM Higgs massj. We fix some of the 
parameters, A , /i, tan/3, and m A to reduce the number of parameters, m A 
and tan/3 are important to realize small ZZh coupling, because CP-even 
Higgs mass matrix is roughly given by 



m A — (m A + m 2 z ) sin (3 cos (3 

— (777/4 + m l) sin /3 cos /3 m 2 z + 8m 2 H 



(4) 



when tan (3 ^> 1. Here, 5m 2 H is the loop correction to m 2 H which can be large 
due to the large top Yukawa coupling. When m 2 A < m 2 z + 8w? H , the main 
component of the lightest Higgs h becomes H^ which has only very small 
ZZHd coupling. Moreover, when m z + Sm 2 Hu — m 2 A ^ (777/4 + m l) s i n cos A 
h includes only small component of H u . On the other hand, if the off- 
diagonal element is large (i.e., tan/3 is small) and/or the difference of the 
diagonal element is small (i.e. m^ is large or the loop correction is small), 
the H u component in h becomes large, which results in large ZZh coupling. 
Therefore, as discussed in [B], 7 < tan/3 and 90 GeV< m^ < 110 GeV are 



1 If we adopt the nonuniversal sfermion masses, then the naturalness requires that only 
the masses of the stops must be around the weak scale and the other sfermion masses can 
be taken much larger than the weak scale. Eq GUT with horizontal symmetry naturally 
obtains such nonuniversal sfermion masses, and our discussion in this paper can be applied 
to the nonuniversal sfermion mass model [11] , 



required (If CP-even Higgs mass is smaller than 90 GeV, the Z — > Ah process 
gives a severe constraint by LEP experiments.). In this paper, we take rather 
large tan (3 and small tua as tan/3 = 15 and iiia = 96 GeV which satisfy 
that the ZZh coupling is smaller than a half of the SM ZZh coupling in the 
parameter region discussed in this papero We examine two cases for A , 
A = and 250 GeV. For the \x parameter, we have a strong reason to take 
it as the weak scale in the scenario with the small ZZh coupling. Since the 
both mass parameters m 2 = m 2 H + |yu| 2 and m\ = m 2 Hu + |//| 2 in the Higgs 
potential must be around the weak scale to obtain the small ZZh coupling, 
not only the tuning for m\ but also that for m 2 are required if the \x parameter 
is much larger than the weak scale. Therefore, we take fi = 275, 300, 325, 
350 GeV, and for comparison, we examine the case \i = 600 GeV, which is 
not so naturajj. 

To obtain the low-energy parameters (sfermion masses, A term and gaug- 
ino masses) from the GUT scale parameters, we use the 1-loop renormaliza- 
tion group equations (RGEs) from the GUT scale down to the electroweak 
scale. In the calculation, we choose the GUT scale Higgs masses in or- 
der to realize m^ = 96 GeV, tan/? = 15, \i = 275, 300, 325, 350 GeV and 
niz = 91.18 GeV at the electroweak scale. Concretely speaking, we choose 
the GUT scale parameters by calculating the RGEs iteratively to realize the 
parameters we fixed at the low energy. Then we calculate the neutralino 
relic density by using micrOMEGAsl.3.7 [TBI [T4"] package by inputting the 
parameters at the electroweak scale. Here we took the on-shell top mass 
M t = 172.6 GeV [15]. 

We display in FigflJ-FigJllthe relic density of the lightest neutralino in the 
natural light Higgs scenario. In the figures, the light gray area is the cosmo- 
logically preferred region where the neutralino relic density is consistent with 
eq.([2]). The regions with larger and smaller relic density are painted black 
and white, respectively. The horizontal-striped region is excluded because 
the LSP becomes stau, which is a charged particle. The dark gray areas are 
allowed regions for b — ► 57 constraint. Since the observed branching ratio for 
the process, Br{b — > sj) exp , = (355 ± 26) x 10~ 6 [IB] , is now in agreement 
with the SM estimations, Br(b -* sj) 8M = (315 ± 23) x 10" 6 [17J, Br(b -»• 
s 7 ) SM = (357 ± 49) x 10" 6 [IE], and Br{b -> s~i) SM = (298 ± 26) x 10~ 6 
[19j . we seek the region where the MSSM contributions for the process are 
moderate. Here we assume the Minimal Flavor Violation (MFV) and the 
primary contributions coming from the SM, charged Higgs and chargino are 



2 We checked by FeynHiggs2.6.4 [U] that the ZZh coupling is smaller than half of the 
SM coupling in the parameter region we took in this paper. 

3 We calculated for /1 = 250 GeV and found that the neutralino relic density is smaller 
than the observed value in all natural parameter space. 
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Figure 1: A = GeV, /i = 275 GeV Figure 2: A = GeV, [i = 300 GeV 
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Figure 3: A = 250 GeV, \i = 275 GeV Figure 4: A = 250 GeV, /i = 300 GeV 

The (rai/2, mo) planes for A = OGeV (upper), A = 250 GeV (lower), 
fi = 275 GeV (left), // = 300 GeV (right) with tan (3 = 15 and m A = 96 GeV. 
Each light gray area is the region where the relic density is consistent with 
the current observation. The relic density is larger (smaller) than that of 
the light gray area in each black (white) region. In the area with horizontal 
stripes, stau is the LSP. Each dark gray area is the region where Br(b — > sj) is 
consistent with the experiment in the sense described in the text. Dashed line 
is the contour on which the mass of the heaviest CP even Higgs is 114.4 GeV. 



taken into account using input parameters and RG method described in the 
previous paragraph. For simplicity, we require the effective Wilson coefficient 
C-j at b quark mass scale to be within 20 percents difference from the SM 
prediction in the leading order approximation for the process. That is to 
say, the coefficients evaluated at the electroweak scale from 1-loop diagrams 
are translated into that of b quark scale (/if, = 4.7 GeV) values using 8x8 
evolution matrix calculated at 2-loop level [HI [20] • Dashed line denotes the 
heavy Higgs mass bound, m# = 114.4 GeV. In the model with small ZZh 
coupling, the ZZH coupling becomes almost the same as the SM value, and 
therefore, the LEP constraints to the SM Higgs mass can be roughly applied 
to the heaviest Higgs mass in the MSSM. In all parameter region in these 
figures, the lightest CP-even Higgs mass is larger than 90 GeV. If tan /3 is 
fairly large, the constraint from B s —>■ fi + n~ process must be taken into ac- 
count [211 l22l [23] . However, the constraint can be negligible in the parameter 
region we took in our calculation, when tan/? ~ 15. 

There are two reasons for the small relic density. For roughly va\ii < 200 GeV, 
the cross section of processes XiXi -^ A ^ bb and XiXi — ► Z — >■ 66 become 
so large that the neutralino relic density is smaller than the observed relic 
density of dark matter. The neutralino is roughly half as heavy as the CP- 
odd Higgs and the Z boson, so the sum of the masses of two neutralinos is 
nearly on the pole of the CP-odd Higgs and the Z boson in this region. The 
cross section of this process decreases as the gaugino mass grows up, because 
the sum of the two neutralinos masses becomes away from the poles. In the 
CMSSM, this left preferred areas are mostly excluded by the LEP bound. 
On the contrary, all preferred regions are allowed by the Higgs mass bound 
in the models with the small ZZh coupling. There is another area in which 
the neutralino relic density becomes smaller than the observed value when 
the gaugino mass becomes larger. This is because the modes of x?x? — > 
two bosons such as x?X? ~~ * hA, x?Xi — * HA, and x?x? ~^> W ± H T modes 
open as well as x?Xi — y Zh. It is essential that in the natural light Higgs 
scenario, all Higgs bosons are light. In the CMSSM, the relic density of the 
region corresponding to this region is larger than cosmologically preferred 
range [21]. We can plot similar graphs even in the models with small ZZh 
coupling in the case of large /i parameter as we show in FiglS] and FigJU] in 
which /i = 600 GeV. This is because the sum of the cross sections of the 
processes x?Xi ~~ > two bosons becomes too small to obtain sufficiently small 
relic density of the neutralino, because the Higgsino components of the light- 
est neutralino become smaller. In the CMSSM, in addition to the difficulty 
in realizing small /i, the modes XiX? ~^ two bosons except for x5x? ~* Zh 
does not open in small mi/2 region because of the heavier Higgs sector. In 
the black region, the main mode is x^x? — ► bb. We comment on the mo 
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dependence of the relic density. When the sfermion mass becomes smaller, 
XiX? — > two leptons process becomes larger because the cross section of the 
process via slepton t-channel exchange increases. Therefore, the smaller mo 
leads to the larger annihilation cross section and the smaller energy density 
of the lightest neutralino in the region. 

The graph of the relic density depends on the \x parameter strongly as we 
commented. The distance between two allowed band region becomes wider 
in (j, — 300 GeV than in /x = 275 GeV. When /x is larger, the left preferred 
band moves to the left because the lightest neutralino becomes heavier. The 
right preferred band moves to the right, because the coannihilation cross 
sections to two bosons becomes smaller for the larger /x parameter. In this 
scenario, the relic density of the neutralino does not depend on the gaugino 
mass so much when the gaugino mass is larger than 300 GeV in the parameter 
region we scanned. Usually when the gaugino mass increases, the total cross 
section decreases. However, in our scenario, the Higgsino components of the 
neutralino increase, and thus the two effects can almost compensate each 
other. In FigJT] and FigJH] we can explicitly see the mild change of the relic 
density as the fairly broad prefferrered region when /x = 325 and 350 GeV. 

The relic density does not change so much if we enlarge Aq. On the 
other hand, the cross section of the b — > 57 process depends on Aq as in 
the FigfTJ-FigSl There are reasonable regions which are consistent with the 
observed relic density of dark matter and experimental constraints when 
A = 250 GeV. When A Q = GeV, there are no or absolutely thin preferred 
region. However, we do not take this allowed region for the b — > sj constraint 
seriously. This is because the allowed region can be changed if there is other 
contributions to b — > sj as sizable gluino contribution. Moreover, the exper- 
imental value of the Br(b — > sj) is larger than the SM prediction. In order 
to increase the Br{b — > sj), larger SUSY breaking scale is required because 
the chargino contribution decreases the branching ratio. This requirement 
makes the allowed region move to upper right. 

3 Conclusions 

We have studied the thermal relic density of the neutralino in the MSSM with 
the light Higgs sector and reasonable SUSY breaking parameters. Actually, 
we took all the dimensionful parameters as order the weak scale. 

The neutralino relic density with the light Higgs bosons is totally differ- 
ent from the well-known result in the CMSSM. In the natural light Higgs 
scenario, the neutralinos to two bosons processes, such as XiXi - > hA and 
X1X1 ~~ > HA open even when the gaugino mass is small. Furthermore, they 
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Figure 5: A = GeV, fi = 600 GeV Figure 6: A = 250 GeV, fj, = 600 GeV 

The (rai/2, m ) planes for A = OGeV (left) and A = 250 GeV (right) with 
fj, = 600 GeV, tan/5 = 15 and m^ = 96 GeV. The usage of each color and 
line is the same as the previous figures. In these figures, the constraints from 
b — > sj are not presented. 
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Figure 7: A = 250 GeV, fi = 325 GeV Fi S ure 8: A o = 250 GeV, \x = 350 GeV 

The (mi/2, m ) planes for /i = 325 GeV (left) and n = 350 GeV (right) with 
A = 250 GeV, tan/3 = 15 and iua = 96 GeV. The usage of each color and 
line is the same as the previous figures. 



dominate the total cross section because the \x parameter is so small that the 
Higgsino components of the lightest neutralino become comparatively large. 
For those reasons, the relic density become smaller than the observed value. 
In contrast, it is mostly larger than the observed value in the CMSSM be- 
cause all Higgs bosons except for the SM-like Higgs are heavy. The region 
with small relic density cannot be excluded because it is not inconsistent 
if there is other dark matter sources. Furthermore, there are cosmologically 
preferred regions in this scenario, which are not excluded by the experiments. 
Thus the MSSM with the natural light Higgs sector is a good model not only 
for naturalness but also for cosmology. 

The cosmologically preferred regions which we studied so far can be tested 
by future direct searches for the weakly interacting massive particles (WIMP) 
because the small \x parameter makes the spin-independent interaction be- 
tween Higgs and the lightest neutralinos large J25J [26] , unless recent direct 
searches for the WIMP such as CDMSE [27] and XENON10 [28] have ex- 
cluded the regions. In the region where the gaugino mass is roughly larger 
than 300 GeV, even if the relic density is smaller than the observed relic den- 
sity of dark matter, it is not so small. (It is larger than 20% of the observed 
value in FigJTJ-FigJBl) Therefore, even if the main component for dark matter 
around the galaxies has only super weak interaction and cannot be found 
in the direct searches of the WIMP, the searches can find the signal for the 
neutralino which is subdominant component, though the concrete prediction 
becomes difficult. We think it an interesting future subject to study the 
direct detection of the WIMP in our scenario. 
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